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Abstract 

We study the structure of Yukawa couplings in intersecting D6-branes wrapping a fac- 
torizable 6-torus compact space T 6 . Models with MSSM-like spectrum are analyzed and 
found to fail in predicting the quark mass spectrum because of the way in which the 
family structure for the left-handed, right-handed quarks and, eventually, the Higgses is 
'factorized' among the different tori. In order to circumvent this, we present a model 
with three supersymmetric Higgs doublets which satisfies the anomaly cancellation con- 
dition in a more natural way than the previous models, where quarks were not treated 
universally regarding their branes assignments, or some particular branes were singled out 
being invariant under orientifold projection. In our model, the family structures of all 
the standard model particles arise in one of the tori and can naturally lead to universal 
strength Yukawa couplings which accommodate the quark mass hierarchy and the mixing 
angles. 



1 Introduction 

Uncovering the nature and origin of the fermion families and the observed pattern of 
fermion mass hierarchies and mixings is one of the most fundamental issues in high energy 
physics. In the framework of the standard model (SM), the vacuum expectation value 
(VEV) of the Higgs field responsible for the electroweak symmetry breaking generate the 
fermion masses through Yukawa couplings. However, the SM does not address the origin 
of these couplings, and the observed values for the fermion masses are considered as initial 
'input' parameters [1]. In addition, the electroweak Cabbibo-Kobayashi-Maskawa (CKM) 
mixing matrix arising from the matrices that diagonalize the up- and down-quarks mass 
matrices is determined experimentally to have, again, a hierarchical structure [1] where 
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the third generation mixing is mostly with the second generation rather than the first. 
Something similar of hierarchies and mixings happens for the neutrinos and huge efforts 
were done in order to understand this 'flavor problem' of the structure of the fermion 
masses and mixing. Phenomenological studies considered "textures" [2] in the form of the 
mass matrices leading to approximately correct relations, and attempts to understand the 
presence of such "textures" then followed in different flavor models [3] or within Grand 
Unified Theories (see [4] and references therein). 

Despite the insight which can be gained from these phenomenological studies of the 
fermion mass matrices, arguably the true resolution to the flavor problem lies in the do- 
main of the underlying fundamental theory of which the SM would be the low energy 
effective theory. Since at present Superstrings/"M" -theory is the only candidate for a 
truly fundamental quantum theory of all interactions, studies of the flavor structure of 
the Yukawa couplings within four-dimensional superstring models are well motivated. In 
particular, the couplings of the effective Lagrangian in superstring theory are in prin- 
ciple calculable and not input parameters, which allows to address the flavor problem 
quantitatively without introducing ad hoc assumptions. Indeed, the structure of fermion 
masses has been studied in a number of semi-realistic heterotic string models such as 
abelian Z n orbifolds [5, 6] which have a beautiful geometric mechanism to generate a 
mass hierarchy [7, 8] and the resulting renormalizable Yukawa couplings can be explicitly 
computed [9, 10] as functions of the geometrical moduli. An important result of such 
studies was to demonstrate that the trilinear superpotential couplings at the string scale 
are generally either zero or 0(1), such that they can provide a natural explanation for 
the top quark Yukawa coupling [11], while other mechanisms utilizing higher-dimensional 
non-renormalizable operators do generate the lighter Yukawa couplings [11, 12]. 

With the advent of Dirichlet D-branes in type II and type I string theory, the phe- 
nomenological possibilities of string theory have widened in several respects. Type I and 
Type IIB orientifold models [13, 14, 15, 16], where the gauge groups of the effective low 
energy Lagrangian arises from sets of coincident D branes and where the matter fields 
arise from open strings which must start and end upon D-branes, were proposed and 
investigations into their general phenomenological features have been possible. In [16, 17] 
a classification of the matter fields has been extracted based on general grounds and for- 
mulas for the soft terms and renormalizable Yukawa couplings were derived. This has 
enabled a number of studies for the patterns of soft breaking parameters [16, 18, 19] and 
Yukawa textures [20, 21]. Nevertheless, the study of the structure of the renormalizable 
Yukawa matrices and its viability within these scenarios of D-branes at singularities has 
proved to be unable to explain the experimental data, since they would generally lead to 
a variant of the "democratic" texture of Yukawa, and one has to break this democracy 
by perturbative higher order effects or non-renormalizable operators, the nature of which 
are still unclear [21]. However, recent studies of the flavor problem within 'interscting D- 
branes' models [22, 23, 24, 25, 26] seemed more promising [27, 28]. In these models, chiral 
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fields to be identified with SM fermions live at different brane intersections and there is 
a natural origin for the replication of quark- lepton generations. In fact, most models are 
toroidal or orbifold (orientifold) compactifications of Type II string theory with Dp-brane 
wrapping intersecting cycles on the compact space, and typically the branes would inter- 
sect a multiple number of times giving rise to the family structure. Moreover, Yukawa 
couplings between three chiral fields arise from open string instantons stretching a world- 
sheet with triangle shape in whose vertices lie the chiral fields. Each worldsheet contribute 
semiclassically to the Yukawa coupling weighted by e~ A , where A is the worldsheet area. 
This exponential weighting makes very natural the appearance of hierarchies in Yukawa 
couplings of different fermions with a pattern controlled by the size of the triangles. 

Yet, the simple model presented in [27] which is based on D6- branes wrapping cycles on 
an orientifold of T 2 x T 2 x T 2 and has the chiral spectrum of the minimal supersymmetric 
standard model (MSSM) does not really give acceptable fermion masses. It leads to a 
mass spectrum of two massless and one massive eigenvalue for the Yukawa matrices. This 
reproduces the leading effect of one generation being much heavier than the other two 
and, thus, should be considered only as a starting point for a deeper phenomenological 
description where small variations on the set up might give rise to smaller but non- 
vanishing masses for the rest of quarks. In fact, we trace this problem of mass degeneracy 
to the factorizable form T 2 x T 2 x T 2 of the compactified space and to the fact that the 
family structures of the quark doublets, the quark singlets and the Higgses arise, each, 
in one of the tori different from the others. This leads to a Yukawa matrix of special 
'factorizable' form (a,ibj) which has always two vanishing eigenvalues. 

We argue in this paper that one can get more interesting Yukawa structures assuming 
three generations of supersymmetric Higgses {Hf, Hf) i=lj2 ,z- This allows the generation 
of family structures for the quark doublets, singlets and Higgses to take place in only 
one of the tori. In fact, several models with three families, including Higgses, have been 
constructed [30, 31] and were favored by the unification of the gauge couplings in heterotic 
string. Importantly for our analysis, having three families of Higgses would allow easily to 
satisfy the RR tadpole cancellation condition which requires the number of fundamentals 
to equal that of antifundamentals for SU(2). This is because the Higgses can account 
for the six SU(2) doublets needed to be added to the three lepton left-handed doublets 
in order to equal the nine antidoublets of the three families of chiral left-handed quark 
color triplets 3(3, 2). This offers a natural solution to the anomaly cancellation condition 
without the need to put the left-handed quarks in different brane intersections [24, 29], or 
to assume some specific properties satisfied by some of the branes [27]. Moreover, having 
three Higgs doublets introduces more Yukawa couplings which introduces more flexibility 
in the computation of the mass matrices [31], hence one can accommodate the observed 
quark masses and their mixing angles. 

The structure of the paper is as follows. In the next section we briefly review the 
different models for intersecting branes leading to an MSSM-like spectrum and state what 
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gauge symmetry they have. Following this, we describe our model of three supersymmetric 
Higgses and the way it satisfies the anomaly cancellation condition. A brief discussion on 
how to determine the string scale in this class of models is presented in section 3. Section 
4 is devoted for a detailed analysis for the quark masses and mixings. Our conclusions 
are given in section 5. 

2 Intersecting Branes Models 

In this section we start with reviewing the construction of MSSM-like models from in- 
tersecting D-branes. We also set some notation that we will use through the paper. 
The intersecting D-branes scenario offers an interesting way to get chiral fermions. Con- 
sider a bunch of N Dp-branes and another set of M Dp-branes (p > 3) both containing the 
Minkowski space and intersecting at some angle in the (p — 3) extra dimensions. One then 
gets massless chiral fermions transforming as (N, M) under the gauge group U(N) x U (M) 
which allows to represent the SM fermions. In addition, if the extra 6 dimensions are com- 
pact, the intersection of a couple of branes is in general multiple and the replication of 
generations is natural. Recently, a particularly interesting class of models yielding 'just' 
the massless fermion spectrum of the SM was constructed [24, 27] . These models consider 
D6-branes in type IIA string theory compactified on a factorizable 6-torus T 2 x T 2 x T 2 . 
One can wrap a D6-brane on a 1-cycle of each T 2 so it expands a three-dimensional cycle 
on the whole T 6 . We denote the wrapping numbers of the D6 a -brane on the i-th T 2 by 
(n l a ,m l a ). If one minimizes the volume of these three-cycles in their homology class, they 
are described by hyperplanes quotiened by a torus lattice and this implies that the num- 
ber of times two branes D6 a and D6b intersect in T 6 is given by the signed intersection 
number 

hb = {n\m\ - m\n\)(n 2 a ml - m 2 a n 2 b ){n 3 a m 3 b - m 3 a n 3 b ) (1) 

In addition to this, one performs an 'orientifold' projection on the torus represented 
by the product fixi?, where Q is the worldsheet parity operator and R is the reflection 
operator with respect to one of the axis of the tori. The set of fixed points under Q x R 
forms an orientifold plane, namely a subspace of spacetime where the orientation of the 
string can flip. This set has 8 components and corresponds to 06-planes wrapped on 
the 3-cycle with wrapping numbers (nj,mj) = (1,0). Now each D-brane a has a mirror 
image under Q x R denoted by a*. If the brane wraps a cycle [ILJ = (n l a ,m l a ) i= i t 2,3 and 
€ (i) re p resen ts the transversal distance of the brane a from the origin in the i-th torus 
in clockwise sense from the direction defined by [n a ], then the mirror image brane a* 
would wrap a cycle [H a *] = (rfa, ~ m a)i=i,2,3 for rectangular tori, and the corresponding 
translation shift from the origin in the i-th torus is given by e^l = — . 

There were, so far, two ways of embedding the standard model gauge group into 
products of unitary and symplectic gauge groups, and both ways used four stacks a, b, c, d 
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(and their orientifold mirrors) of D6-branes called respectively the Baryonic, Left, Right 
and Leptonic branes. Both methods (models) succeed in getting an MSSM-like spectrum 
free of anomalies. However, in order to do so, in the first model two left-handed quarks 
were doublets and one left-handed quark was an antidoublet, while in the second model 
one of the branes (b) was singled out being invariant under the orientifold projection. We 
will summarize in the next subsection the set up of these two models then we will present 
in the following subsection our set up to generate the SM-like spectrum with the aid of 3 
supersymmetric Higgs doublets. 

2.1 Models with MSSM-like spectrum 

In the first model (see [24] for details), and let's call it model A, one gets initially the 
gauge symmetry 

Model A : 2/(3) x 17(2) x 17(1) x 17(1) (2) 

resulting from the following number of branes in the corresponding stacks: 

N a = 3,N b = 2,N c = l,N d = l 

and then one would embed SU(3) C into 2/(3) and SU(2) L into U(2). In order to yield the 
desired SM spectrum, it is enough to select the wrapping numbers (n l a , m l a ) for the four 
sets of D6-branes in such a way that the intersection wrapping numbers are given by: 

hb = 1 i lab* = 2 

T — — 3 ■ T — —3 

A ac ° : A ac* u 

hd = ; hd* = —3 

hd = —3 ; hd* = 3 (3) 

all other intersections vanishing. The massless fermion spectrum living at the intersections 
is shown in Table 1, where the N R represents a right-handed neutrino and the hypercharge 
generator is defined as: 

Qy = \Qa ~ \Qc + \Qd- (4) 



In this model one adopts the choice of splitting the left-handed quarks into one quark 
represented by the intersection (ab) and the other two reprsented by (ab*) for consistency 
requirements. In fact, as already mentioned, the RR tadpole cancellation condition, which 
is stronger than the gauge anomaly cancellation condition, requires the same number of 
doublets and antidoublets. This choice, then, allows the left-handed quarks not to be 
universal under the 27 (1)& charge so that if two left quarks were 2/(2) doublets and the 
other one was U(2) antidoublet, then taking the SM leptons as U(2) antidoublets allows 
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Table 1: Standard model spectrum and U(l) charges in the first model (A) 

to satisfy the requirement without the need of extra doublets. As to the Higgs field sector, 
the Higgses would come from the intersection between b(b*) and c(c*) branes, and there 
are four possible varieties of them (hi, Hi)^^ since we have two varieties of left quarks 
(Ql,Ql) and two varieties of right quarks (Ur,Dr). 

The second model (see [27, 28] for details), to be called model B, presents a slight 
variation where N^ — l but b* , the mirror of b, lies on top of it (b = b*), so it can actually 
be considered as a stack of two branes which, under Q projection, yields a U Sp(2) = SU (2) 
gauge group. So the initial gauge group is 

Model B : 17(3) x SU(2) x U(l) x U(l) (5) 

With the following intersection numbers 

lab = 3, 
lac = 3, 
Idb = 3, 
Idc = —3, 

and all the other intersection numbers being zero, one gets the spectrum shown in Table 
2. The Nr denotes the right-handed neutrino and the hypercharge generator is defined as 
Qy — \Qa — \Qc — \Qd- Notice that we do not have here the Qb anomaly condition since 
doublets and antidoublets in SU (2) are the same (there is no U(l\ to differentiate between 
them). A particular class of configurations satisfying the conditions (6) is presented in 
Table 3 and the intersections (be) and (be*) can be identified with the MSSM Higgs 
particles H u , H d . 

2.2 Models with 3 Supersymmetric Higgs doublets 

As discussed above, the first model (A) treats the left quarks differently as regards their 
location at the branes intersections. Moreover, the intersection numbers (eq.3) are not 
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Table 2: Standard model spectrum and U(l) charges in the second model (B). 
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Table 3: D6-brane wrapping numbers giving rise to the chiral spectrum of the MSSM in the 
second model (B). 

'symmetric' among the branes and their mirrors (e.g., hd = 0, ha* = 3). As to the second 
model (B), although it also reproduces an SM-like spectrum (with a right-handed neu- 
trino), it singles out one of the branes by requiring its invariance under orientifold action 
(b = b*). The origin of these assumptions lie, as already said, in the consistency require- 
ment that the number of fundamentals should be equal to the number of antifundamentals 
even for SU (2). 

We are proposing now another way to satisfy this condition. We consider, as in the 
first model (A), a stack of two branes (b) giving rise to U(2) gauge symmetry. We will treat 
the three left-handed quarks universally and consider that we have chiral quarks in 3(3, 2) 
under SU(3) x SU(2). The full model must contain then nine fields (1, 2), three of which 
correspond to left-handed leptons. As to the remaining six doublets, we do not need extra 
doublets to be accounted for if we take the natural assumption of three generations of 
Higgses (H™, Hf) i= i2$- In fact, in both models (A) and (B) the u-Higgs and the d-Higgs 
are assigned opposite U(l) charges so that anomaly would not be affected by including 
them. However, no reason prohibits them from having the same U(l)b charge, and thus 
they can provide the extra six doublets necessary for anomaly cancellation 

Our model will be purely toroidal, with no orientifold projection and so no mirror 
branes added. We shall consider that the SM particles reside at the intersections amongst 
4 stacks of branes N a = 3, N b — 2, N c — 1, N d — 1. The intersection numbers would be 
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then 



\I a b\ = 3 representing Ql , \hc\ = 3 representing H u , 

\I ac \ = 3 representing Ur , \Ibd\ = 6 = 3 + 3 representing H d ,L , (7) 

\I ad \ = 3 representing , |7 cd | = 3 representing E R , 

In this third model, and let's call it model C, we used the fact that the Higgs H d and 
the lepton L have the same SU(3) x SU(2) x U(\)y quantum numbers, and so one can 
consider getting both of them at the intersection of the same branes (b and d). Requiring 
that the observed hypercharge generator is a linear combination of the four £7(1) 's one 
finds the following general solution 

Qy = (^ + Pa)Q a + (^ + Pa)Q b + aQ c + -f(l + Pa)Q d (8) 

where f3 2 = 7 2 = 1 defined in Table 4 and a is arbitrary. Notice that with the choice 
a — — |, /3 — +1, 7=1, we get the hypercharge defined in equation 4 
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Table 4: Standard model spectrum and £7(1) charges in the third model (C) 

We will give in section 4 an example of D6-branes wrapping numbers realizing the 
conditions (7). As can be seen from Table 4, all the U(l) gauge groups are anomaly 
free. One should consider also mixed gauge and gravitational anomalies. However, we 
expect that these anomalies are cancelled by a generalized Green-Schwarz mechanism 
and that three combinations of the £7(1) 's would get massive with mass roughly of the 
order of the string scale, while the hypercharge Y combination would stay massless. The 
symmetries whose gauge bosons become massive would disappear as gauge symmetries 
from the low-energy effective field theory, but remain as global symmetries unbroken in 
perturbation theory. In this respect, U(l) a and U(l)d represent respectively the global 
Baryonic and Leptonic number symmetries. However, assigning Qd charges to the Higgs 
H d might lead to a breaking of the Lepton number symmetry when the Higgs acquires 
a VEV. Notice that in this model we do not have a right-handed neutrino as a chiral 
fermion from intersecting branes. Also, once we assume the Higgses H u , H d come in a 
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number of generations equal to that of the SM particles then the gauge anomalies would 
be cancelled automatically. Thus, in this model there is no relation between the number 
of colors and the number of families as it was the case in the previous models (A) and 
(B) [27, 24]. 



3 String scale in the model with three Higgs doublets 

The toroidal models are in general non-supersymmetric, and might have tachyons at 
intersections. However, it is possible to vary the compact radii in order to get rid of the 
tachyons. Also, one can adjust the radii so that there is one massless scalar at any given 
intersection, which means that one gets N=l SUSY at that specific intersection [24]. For 
instance in the MSSM-like model (B), if the ratios of radii in the second and third torus 
are equal then the same N=l SUSY is preserved at all intersections and the model is 
locally N=l supersymmetric, having an MSSM spectrum with a minimal Higgs set [27]. 

For non-supersymmetric models, stabilization of the hierarchy of the weak scale can be 
achieved by lowering the string scale down to few TeV [32, 33], while for supersymmetric 
models there are several arguments in favor of string scales in "intermediate" range Mj 
10 10-14 GeV [20]. Such arguments provide an explanation to the observed experimental 
neutrino masses [35], or a means to attack the hierarchy problem of unified theories in 
supergravity models by getting the gravitino mass around the weak scale m 3 / 2 ~ Myy in a 
natural way without invoking any hierarchically suppressed non-perturbative effect [36]. 
Also, for intermediate string scale scenarios, charge and color breaking constraints on the 
acceptable region of parameter space for soft supersymmetry breaking terms become less 
important [37, 38]. In addition, the observed ultra high energy (10 20 eU) cosmic rays can 
be explained, for intermediate string scale, as products of long lived massive string mode 
decays [39]. 

In order to compute the string scale at which the running coupling constants intersect 
in model (C) with three Higgs doublets, one uses the one-loop running equation for the 
gauge coupling 

1 1 +^ + ^. <«» 



a t (Q) Oi(M z ) 2it M z 2tt M s 

where a, = gf/^n with % = 2,3, Y, and ftj's are the coefficients of the /^-functions, and 
where M z represents the overall non-supersymmetric scale while Ms ~ 500GeU represents 
an overall supersymmetric scale [40]. On the other hand, from eq.(8) we have the following 
relation at the string scale Mj 

1 a 2 , (l + (3a) 2 | (l/2 + Paf | (2/3 + a(3) 2 



ay(Mj) af(M/) af(Mj) a 2 (Mj) a 3 (M 7 

For the SM content with one Higgs doublet, the non-supersymmetric /^-functions are given 
by b^ s = 7, b^ s = 19/6 and by S = -C 2 x 41/6, where C is the normalization constant of 
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the U(l)y hypercharge (C 2 = 3/5 in 577(5) GUT). As for the supersymmetric /3-functions, 
considering three supersymmetric generations of standard particles, two Higgs doublets 
and an arbitrary number of extra particles, we have 

bl = 3-±n 3 , (11) 

b! = ~l-\n 2 , (12) 

1% = -C 2 x (11 + q) , (13) 

where 

ni n 2 n 3 

i=l j=l k=l 

and n-i, n 2 , n 3 are the number of extra SU(3) x SU(2) singlets, SU(2) doublets and 
SU(3) triplets, respectively, with masses close to M s and hypercharges Y t . From eqs.(9) 
and (10), one finds 

ln Mj_ = 1 ( , 1 _ a 2 _ (l+{3a) 2 

n M s ~ [(| + a /3)26f + (I + a p) 2 b$-b Y ]\ ^a Y {M z ) of(Mj) of(Mj) 

We shall use the experimental values [41] M z = 91.187 GeV, a 3 (M z ) = 0.1184, a 2 (M z ) = 
0.0338, ay{Mz) = 0.01016 given in the MS scheme. The fact that we have four extra 
Higgs doublets with respect to the case of the MSSM means that we should take n 2 = 4 
and rii = n 3 = 0. Now, assuming that a = 1/2, (3 = +1 and C 2 = 3/5, one obtains 

1- ^ - 40.56 - 4" - 4» - 1.89 In ^. (16) 

Thus for c^(M 7 ) ~ af(Mj) ~ 0.1, one finds M 7 ~ 10 12 GeV. One could also check 
that the curves of a 2 and a 3 cross at approximately this intermediate scale which is, as 
emphasized above, an attractive possibility. However, we should mention here that in this 
class of toroidal or orientifold models there may exist extra chiral fields that would change 
the gauge couplings and might lead to a lower string scale. Although this possibility is 
indeed crucial in non-supersymmetric models in order to avoid any hierarchy between the 
string scale and the electroweak scale, it is not essential in our model with supersymmetric 
content. The model could still be consistent if these additional fields are decoupled from 
our spectrum or the possible threshold corrections are small [33, 34]. 
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4 Analysis of Fermion Masses and Mixing 



In [27], it was shown that a Yukawa coupling between fields at the intersections of factor- 
izable 3-cycles Il a , Tl b and Il c on a factorizable T 6 is given by 



r=l 



Sir) 
0« 



(17) 



Here, each triplet of intersection (i, j, k) is described by the following multi-indices 



^(i(D,i(2)^(3) )e n a nn 6 , 
j = (j (1) ,j (2) ,j (3) )en c nn a , jW = o,...,|/W| 
fc = (few, fc( 2 >, A;( 3 )) e n 6 n n c , fcW = o, . . . , 



i( r ) = n l/^l-l 

1, 
1, 



(18) 



where (r) is an index indicating the r torus, and I, 



ah 



n 



m':'), where lf£ denotes the intersection number of cycles a and b on the r th torus 



n n 7~ (r) 

ll r =l A ab 



n 



Jh 



and is the total intersection number; <r abc = sign(I ab I bc I ca ), h qu stands for the quantum 
contribution to the instanton amplitude and d is the complex theta function 
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e Tri(5+l) 2 K e 2ivi(8+l)<f> 
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We have 
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l ab 



j(r) f,(r) 



T (r) 
i ca 



j(r) 
1 bc 
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r(r) Ar) Ar) 
1 ab 1 bc lca 



I no k ca I 



ft' 



(20) 
(21) 
(22) 



where g^ 7 *** = 5f.cc?. (i^, /5c" 1 ' anc ^ s<r ' ) = s {^ r \j (j \ k^) G Z is a linear function of the 
integers i^ r \ and k^'\ represents the 'shifts' in the r th torus while the phase 9^ 
accounts for adding Wilson lines around the D-brane wrapping 1-cycles in the torus. 
j(r) re p resen t s the Kahler structure of the r^ torus and so would be proportioanl 
to its area. Once, we have determined the Yukawa couplings, one can compute the 
quark masses and mixings to see whether the model reproduces the observed hierarchical 
structure. 



4.1 Models with one supersymmetric Higgs doublet 

In the first model (A) [24], the case of a minimal set of Higgs fields similar to the MSSM 
was shown to give masses to the top, charm and bottom quarks while the strange, down, 
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and up quarks remained massless. It was argued that, with a double Higgs system case, 
the observed hierarchy of fermion masses would be a consequence of the different values of 
the Higgs fields and the hierarchical values of Yukawa couplings, coming from geometrical 
considerations. 

As to the seconde model (B) [27] and with the wrapping numbers shown in Table 3, 
one gets one Higgs doublet H u (H d ) at the intersection of bc(bc*). Yukawa couplings of 
the form 

YVQIH U UI Y?QiH d D R \ (23) 

were computed, and only the third generation of the quarks are massive. It was argued 
that smaller perturbation of this set up can give rise to smaller but non- vanishing masses 
for the rest of the quarks. 

However, examining the model in depth would trace the problem of having two zero 
eigenvalues in the Yukawa matrices to the 'factorizable' form that they take when the 
family replications for the left-handed quarks and the right-handed quarks come from 
different tori. For the case of Table 3 we see that the index = 0, — 1, 



denoting the left quarks, would span the values 0, 1, 2 only in the 2 nd torus, while the index 



l ac(c*) 



= 0,..., 



]1 ac(c*)\ 



1, denoting the right quarks, would have its family structure only 



in the 3 rd torus. In such cases, and neglecting Wilson line effects, the Yukawa couplings 
would always be of the form 



0<D 







(^)x^ 5 f (^)x^ 5 ^ (V 3 )) 







(24) 



and so it is of a 'factorizable' form 



Yij ~ aibj 



Such matrices always have two zero eigenvalues since, for instance, the second and third 
columns are proportioanl to the first one. 

Could we get more interesting phenomenology if the family structures of both left- 
handed and right-handed quarks arise in the same torus? The answer is no, if we restrict 
ourselves to one Higgs doublet. In fact, if we adopt the wrapping numbers shown in Table 
5, where the branes b, c, c* are on top of each other in the second torus, we find that the 
conditions (6) are satisfied, and we have one massless non-chiral Higgs doublet arising at 
the intersection of the brane b and the brane c (or c*) in the first and third tori. In other 
words, there is a minimal Higgs sector with a /z-parameter determined by the distance 
between the branes b and c along the second torus. 

In order to compute the Yukawa structure, we now note that the family structure for 
both the left-handed and the right-handed quarks is originated in the second torus, and 
so, neglecting Wilson lines effect, we shall get 



^(i) 



5(0) 




o 



5(0) 




(«< s >) (25) 
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where 5^ 2 \i, j) = ^ + A, and A is a constant determined by the shifts e a , e^, e c . However, 
using the periodicity of theta function 



0(r) 



5+1 



« 



(26) 



we get the following form for the Yukawa matrix 



a 6 c \ 
6 c a 
c a b j 



(27) 



A matrix of this form has a spectrum such that two of the eigenvalues are always opposite 
in sign, so it leads to two degenerate states and can not reproduce the hierarchy in the 
masses of the quarks. Having spotted the origin of the problem, we now move to our third 
model (C). 



4.2 Model with 3 supersymmetric Higgs doublets 

We saw in the previous subsection that having the family structures of the left-handed 
quarks and the right-handed quarks to arise from different tori leads to a mass matrix 
with two vanishing eigenvalues. Also, having one Higgs doublet in the set up would lead 
to a phenomenologically unacceptable form for the mass matrices. One could also check 
that getting Higgs doublets replication in one torus different from the torus where the 
family structure arises for the quarks, is similar to the one Higgs doublet situation. In this 
case, their effects are factored out. Thus, we are led naturally to seek a situation where we 
have more than one Higgs doublet and where the family structures of all the left-handed 
quarks, right-handed quarks and the Higgses arise in one of the tori. Recalling that the 
assumption of three supersymmetric Higgses is also a 'normal' choice in order to cancel 
anomalies, we consider the model (C) with 3 Higgs doublets. Since we are interested in 
this paper only in the quark sector, we adopt the intersection numbers (7) of the model 





«, ml) 






N a = 3 


(1,0) 


(1,3) 


(1,0) 


N b = l 


(0,1) 


(1,0) 


(o,-i) 


N c = l 


(1,1) 


(1,0) 


(1,1) 


N d = l 


(1,0) 


(1,3) 


(1,0) 



Table 5: Alternative example of D6-brane wrapping numbers in the second model (B) leading to 
a chiral spectrum of the MSSM. The family structure of both the left-handed and right-handed 
quarks arises in the second torus. 
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(C) seeking to generate family structures for the left and right quarks, as well as for the 
Higgses, in the second torus, say. No constraint is imposed on where the family structure 
for the leptons would arise. This means that we are not trying here to interpret the lepton 
masses and mixing, in particular that the model does not contain right-handed neutrinos, 
nor Majorana neutrinos because Lepton number is a symmetry, and so the question of 
neutrino masses should be addressed differently. 





(nl, ml) 


(n 2 a , ml) 


(n 3 a ,m 3 a ) 


N a = 3 


(1,0) 


(2,3) 


(1,0) 


N b = 2 


(0,1) 


(1,3) 


(0,1) 


N c = l 


(0,1) 


(1,0) 


(0,1) 


N d = l 


(1,-1) 


(3,3) 


(1,-1) 



Table 6: Example of D6-brane wrapping numbers in the third model (C). The family structure 
of the standard model particles arises in the second torus. 

The conditions (7) are obtained with the wrapping numbers shown in Table 6 where 
the branes b and c are on top of each other in the first and third tori, and so one gets 
massless non-chiral H u Higgs doublets. Moreover, we see that the family structures of 
all the standard model particles (the left and right handed quarks, the left and right 
handed leptons and the Higgses) arise in the second torus. The branes b and d give 
an intersection number equal to 6, so we can identify the first three intersections as the 
three H d Higgs doublets while the last three intersections would be the three left lepton 
doublets. In this way the Higgs doublets H d , like the left leptons, are chiral but this 
does not lead to a problem in constructing the MSSM superpotential since the chiral 
H d can still form a /i-term with the non-chiral H u . An approach to deal with chirality 
issues is to compactify over T 2 / Z 2 instead of T 2 [23]. However, one should compute the 
detailed spectrum to check when this would generically project onto chiral matter. We do 
not follow this approach here, but instead, seek an assignment of wrapping numbers that 
leads to chiral fermions. We found that such an assignment could be obtained provided we 
allow multiwrapping cycles for the branes. As an example, in Table 6 we use a cycle (3, 3) 
for the brane d in the second torus. Since the wrapping numbers are not coprime the brane 
d is multiwrapped 3 times over the cycle (1, 1) in this torus. Normally, multiwrapping 
leads to an enhancement of the gauge symmetry [42] (look also at [23, 43] where similar 
multiwrapped assignments were used for D4-branes, and hence for the whole compact 
dimension of the brane). Nonetheless, in our case the multiwrapping occurs only in the 
second torus. Even if this partial multiwrapping in the second torus enhanced the world- 
volume gauge group from U(l)d to U(l) d (with generators Q%, a = 1, 2, 3), our discussion 
regarding the hypercharge and the anomaly cancellation (equation 8) would stay valid 
with Q d = J2 a =i Qd [43] • Thus, we shall not examine further the effects of multiwrapping, 
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especially that our model should be considered as a step towards building a more realistic 
one. Actually, the wrapping numbers in Table 6, as it is the case in Table 3 and Table 
5 (see [27, 44] for discussion on this point), show that the corresponding brane content 
by itself does not satisfy the RR tadpole conditions I] a iV a n a = 0, which would read as 
follows 



E a N a n\nlnl = 
EaNamlnln^O 
E a N a nlm 2 a nl = 
J2aN a nln 2 a ml = 



E a N a n\m 2 a ml -- 
Ea N a mln 2 a ml -- 
Ea N a mlm 2 a nl -- 
Ea N a mlm 2 a ml 







= 



(28) 



Yet, since tadpole cancellation conditions are closely connected to cancellation of anoma- 
lies, and our model does cancel the anomalies related to the gauge groups in Table 4, 
it is not so surprising that with a slight change in the set up one could satisfy the RR 
tadpole conditions. In fact, and in the spirit of bottom to top approach, the model should 
be seen as a submodel embedded in a bigger one where extra RR sources are included. 
These may either involve some hidden, and possibly non-factorizable, extra branes with 
no neat intersections with the SM branes or some NS-NS background fluxes, with none 
of these possibilities adding a 'net' chiral matter content [45]. We shall not dwell on 
the details of such embedding which might lead to extra matter, expectedly, heavy and 
disconnected from the SM sector. Rather, we shall take our set up and examine what 
interesting geometrical explanations for the fermion masses it might lead to. The quark 
Yukawa coupling with the H d Higgs would then be proportional to a product of three 
theta functions (neglecting again the Wilson line phase) 



0(U 



5(0) 




(V 1 )) X ^ 



S(i,j, k) 




(29) 



with k = 0, 1, 2. The index k runs only over the first three 1m intersections identified 
with the H d Higgs doublets. Thus the quark Yukawa couplings for both the H u and H d 
Higgses would be proportional to 

8(i,j, k) 




Yijk ~ 



# 2 ) 



(30) 



and so we will restrict, henceforth, our discussion to the second torus. For the U-quark 
Yukawa coupling Y U H U Q L U R we have \I ab \ = \I^\ — \hc\ = 3. This is similar to the case 
of elliptic fibration discussed in [27, 46] where the intersection numbers are not coprime 
and only the triplets of intersection satisfying the selection rule 

% + j + k = mod 3. (31) 

are connected by an instanton. We then get the following Yukawa couplings: 



ijl 




Y 



ij2 



( 








C 







A 





V 


B 
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Y 



ij3 



( o 

c 

\ 



B 





o \ 





(32) 



with 



A = ■& 



e/3 




(3J/c/), B = -d 



(e-l)/3 




(3J/c/), C = tf 



(c + l)/3 




(3J/a') 



(33) 

and where we have e = e a + e& + e c . For the D-quark Yukawa coupling Y d H d QiDR one 
would get the same result with a different e-shift e' = e a + £5 + e d . However, as we shall 
see, a numerically good fit is obtained around e ~ e' ~ 0, and to fix the ideas let's take 
e — e'. Thus, we get the U-quark mass matrix 



(34) 





( Av\ 


Bv% 


Cv u 2 


\ 


Ml = h qu 


Cv 3 










{ Bv u 2 


Cv'l 


Avl 


) 



and the D-quark mass matrix 





( Avf 


Bv d 


Cv d \ 




Cvl 


Av$ 


Bvf 




[ Bv d 2 


Cvf 


Av d ) 



(35) 



where h qu includes the quantum fluctuation factor and we expect it to be similar for the 
u- and d- quarks since leading effects would come from QCD loops [28], and vf ,d is the 



VEV for the Higgs H?' d with 



E«) 2 + K") 2 = (174) 2 (Ge\/) 2 

i=l 

The quark masses are obtained by diagonalizing the above mass matrices 

U L M u Ui = d v 
D L M d rf R = d D 

where Ul, Ur, D l , D r are unitary matrices and 

djj = diag(m t ,m c ,m u ) 

do = diag(m b ,m s ,md) 
while the CKM matrix is given by 

CRM = U L D[ 



(36) 



(37) 



(38) 



We have 7 free parameters consisting of the 6 Higgs VEVs with the constraint (36), the 
area of the torus and the shift e. We will not include the unknown overall multiplicative 
factor h qu which is of order 0(1). This set of parameters can be fixed by the quark masses 
and one mixing angle and the model has to predict the remaining two mixing angles in 



16 



the CKM. This might be a non-trivial task since one has to span the whole range of all 
of these free parameters very carefully. Here we will consider some examples and try to 
show that for a particular choice of these free parameters one may obtain a well studied 
Yukawa texture, like for instance the Universal Strength Yukawa (USY) couplings (see 
Ref.[47], and references therein). Let us start with the case of approximately symmetric 
matrices M u , M d i.e, B m C. In this case, the shape of theta function for a fixed area 
argument shows that it is centered symmetrically around e = 0, and so we will span the 
e-parameter around this value. Also, in order to generate the mass spectrum one could 
put the mass matrices in the form 



M^f = h qu Av u / 



/ u.d i u.d u,d I u.d \ 

( V l / V 3 "1 "2^2 l V 3 \ 

u 
3 

u,d j u,d u,d I u,d 



(39) 



u,a I u,a u,a / u,a 

where ct\ = B/A and a-i — CjA. So one could generate the spectrum provided that 

{<, v%} oc {m u , m c , m t } , 

{vf, vi V3} « {m d , m s , m b } . (40) 

and a 1 ,a 2 1. The conditions (40) with the constraint (36) would determine the range 
in which we should vary the VEVs. With such considerations one finds that the following 
choice of parameters: 

~ 63 MeV, vl ~ 0.95 GeV, v% ~ 174 GeV 
vf ~ 8.5 MeV, vi ~ 136 MeV, v% ~ 4.2 GeV 

e ~ 0.002 , area ~ 18.71 . (41) 

gives the following quarks mass spectrum 

djj = {mt = 173.9 GeV, m c = 1.02 GeV, m u = 4.3 MeV} 

d D = {m b = 4.19 GeV, m s = 136 MeV, ra d = 8.2 MeV} (42) 

which are in the experimentally acceptable range [1], and a CKM matrix with diagonal 
elements near the unity, and (Vo km) 12 — 0.216. However, (V C km)i3 ~ (V C km)23 ~ 
10~ 4 - 10- 5 . 

One can also look for other structures different from the "hierarchical" (a <C 1) tex- 
ture, for example the case of cti^ ~ V1/V3 ~ ^2/^3 ~ 1 leads to a nearly democratic 
Yukawa texture which is known to accommodate the observed masses and mixing. How- 
ever, our checks with real-valued vevs indicate that our configuration leads to the correct 
masses and one mixing angle only while the other two mixing angles are smaller than 
their experimental values. This is similar to Z3-heterotic situation in [31], where another 
mechanism, Fayet-Iliopoulos breaking, was called for in order to address the question of 
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the complete quark mixing. However, with complex vevs and a democtratic texture, one 
gets the following USY texture. 




- u.d \ 

1 e^23 \ 

• u,d ■ u.d 

- u.d 

e VPl3 1 / 



(43) 



This type of Yukawa, where all Yukawa couplings have the same modulus and the flavor 
dependence being all contained in the phases, has been recenetly studied in [47]. It was 
shown that with very small values of the phases ~ 1CT 3 — 1CT 2 one could generate the 
right values of the quark masses and the CKM mixing angles. It is interesting to note 
that this class of couplings is motivated by horizontal symmetries [47] and also arises in 
the models with two large extra dimension [48]. Here we find another motivation for the 
USY couplings. 

5 Conclusions 

We have shown how simple configurations of D-branes wrapping a compact space may give 
a good quantitatively description of quark masses and mixing. In particular, one finds that 
with a 3 supersymmetric Higgs doublets model the anomaly cancellation condition could 
be solved easily without introducing extra matter doublet fields, nor putting assumptions 
on the quarks' brane assignment or on the branes themselves. In this class of models, 
it turns out that the string scale is of order 10 12 GeV which is an interesting scale for 
generating neutrino masses and many other phenomenological issues. With real Higgs 
vevs, the model can easily account for the quark masses and one of the CKM mixing 
angles. However, with complex vevs one can get Yukawa couplings in the form of USY 
textures which can accommodate the masses and the three CKM mixing angles with very 
small phases. It would be worthwhile to study the leptonic sector in this perspective. 
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